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SUMMARY 


The  problem  of  the  minimum  induced  drag  of  wings  haying  a  given  lift 
and  a  given  span  is  extended  to  include  cases  in  which  the  bending  moment 
to  be  supported  by  the  wing  is  also  given,  “As  in  the  classical  problem 
of  induced  drag,  the  theory  is  limited  to  lifting  surfaces  traveling  at 
subsonic  speeds.  It  is  found  that  the  required  shape  of  the  downwash  dis¬ 
tribution  can  be  obtained  in  an  elementary  way  which  is  applicable  to  a 
variety  of  such  problems.  Expressions  for  the  minimum  drag  and  the  cor¬ 
responding  spanwise  load  distributions  are  also  given  for  the  case  in 
which  the  lift  and  the  bending  moment  about  the  wing  root  are  fixed  while 
the  span  is  allowed  to  vary.  The  results  show  a  15— percent  reduction  of 
the  induced  drag  with  a  15— percent  increase  in  span  as  compared  with 
results  for  an  elliptically  loaded  wing  having  the  same  total  lift  and 
bending  moment. 


INTRODUCTION 


In  the  problem  of  minimum  induced  drag  as  originally  treated  by  Munk 
(references  1  and  2)  the  span  of  the  wing  and  the  total  lift  were  supposed 
to  be  given  and  the  distribution  of  lift  over  the  span  resulting  in  a  min¬ 
imum  of  drag  was  sought.  The  solution  of  this  problem  thus  provided  a 
convenient  lower  bound  for  the  induced  drag  of  a  wing  of  given  dimensions. 

In  the  practical  design  of  wings  the  requirements  for  low  induced 
drag  and  the  requirements  for  structural  strength  are  opposed.  Here  the 
bending  moment  developed  by  the  lift  becomes  an  important  considera¬ 
tion  —  more  important  in  many  cases  than  the  actual  spanwise  dimension 
of  the  wing.  Such  considerations  lead  to  the  problem  of  determining  the 
minimum  drag  with  limitations  imposed  on  the  bending  moment  as  well  as 
on  the  total  lift.  It  is  the  purpose  of  the  present  paper  to  show  how 
the  methods  of  the  earlier  analysis  can  be  extended  in  a  very  simple  way 
to  the  solution  of  problems  involving  the  bending  moment  of  the  load 
distribution. 


2 


MCA  TN  22^9 


A  complete  list  of  symbols  employed  in  the  analysis  will  he  found  in 
the  appendix. 


GENERAL  FORMOLAS  FOR  LIFT,  DRAG,  AND  BENDING  MOMENT 


Reference  may  he  made  to  the  original  papers  of  Prandtl  and  Munk 
(references  1  and  2),  or  to  any  of  the  standard  text  hooks  on  aerodynamics, 
for  the  fundamental  developments  of  wing  theory  which  form  the  basis  for 


the  calculations  of  induced  drag, 
is  given  hy 


and  the  drag  is  given  hy 


L  =  pV 


Dp  =  p 


In  these  developments  the  over-all  lift 


+s 


T  dy 


— s 


+s 


w,-  r  dy 


(1) 


(2) 


— s 


In  these  formulas  the  wing  span  is  supposed  to  extend  along  the  y  axis 
between  — s  and  +s,  T  is  the  local  circulation  or  vortex  strength,  and 
V  is  the  constant  velocity  of  flight.  The  induced  downwash  velocity  Wj_ 
is  variable  along  the  span  and  is  connected  with  the  vortex  distribution 
f  (y)  through  the  relation 


*i(y)  =  57 


y-*i 


(3) 


With  this  value  for  wi  the  expression  for  the  drag  may  be  converted  to 
a  double  integral  involving  the  spanwise  distribution  of  lift  as  repre¬ 
sented  by  the  circulation  strength  F 


*  ■  hr  l 


+s  p+  s  r(y)r’(T)) 


y-*i 


dy  dq 


w 


This  integral  may  be  reduced  to  a  more  symmetric  form  if  it  is  integrated 
by  parts  on  the  supposition  that  F  falls  to  zero  at  the  wing  tips.  Thus1 


Di  = 


r(y)r(n) 

(y-ri)2 


dy  dr) 


i-The  validity  of  equations  (3),  (4),  and  (5)  can  he  demonstrated  hy 

referring  to  the  limiting  values  of  complex  integrals  taken  along  a  path 
a  short  distance  above  the  singular  point  on  the  real  axis.  In  the  case 
of  equations  (3)  and  (4)  this  process  yields  the  Cauchy  principal  value. 


NACA  TN  22^9 


3 


In  mathematical  terms  the  problem  is  to  minimize  the  double  integral, 
equation  (5)>  while  holding  fixed  values  of 


L 


PV 


p+S 

J-S 


r  d y 


and 


(y-s0)  r  <*y 


(7) 


where  B  is  the  bending  moment  specified  about  the  point  s0.  For  the 
time  being  s0  will  be  taken  as  the  origin,  or  wing  root  (s0  =  0), 
although  later  another  example  will  appear. 

Although  the  derivation  of  the  formulas  for  induced  drag  makes  use  of 
the  concept  of  the  lifting  line,  it  is  important  to  note  that  the  results 
are  not  actually  restricted  to  this  approximation.  According  to  Munk?s 
well-known  stagger  theorem  the  induced  drag  of  a  lifting  surface  will  be 
equal  to  that  of  a  lifting  line  if  the  spanwise  load  distributions  are 
the  same. 

It  should  be  noted  further  that  the  induced  drag  of  a  wing  having  a 
given  lift  and  a  given  spanwise  load  distribution  is  not  affected  by  the 
compressibility  of  the  air  at  subsonic  speeds.  At  supersonic  speeds  an 
additional  drag  associated  with  the  formation  of  waves  arises  and  the 
induced  drag,  which  is  associated  with  the  vortex  wake,  becomes  only  a 
part  of  the  total  pressure  drag. 


THE  DISTRIBUTION  OF  DOWNWASH  FOR  MINIMUM  DRAG 


In  general,  if  the  drag  is  to  be  a  minimum,  a  small  variation  in  the 
shape  of  the  curve  of  spanwise  loading  will  produce  no  first-order  change 
in  the  drag.  The  variation  in  shape  may  take  the  form  of  a  small  addition 
to  the  original  loading;  it  is  then  necessary  to  find  conditions  under  . 
which  the  drag  added  by  a  small  additional  loading  is  zero. 

The  solution  of  this  latter  problem  is  rendered  especially  simple 
by  the  mutual  drag  theorem  (reference  l),  which  arises  from  the  evident 
symmebry  of  the  integral  to  be  minimized  (equation  (5))*  The  theorem 
states  that  if  the  lift  distribution  (represented  by  T  (y))  is  the  sum 
of  two  distributions  T\,  and  T2,  the  drag  of  T x  arising  from  the 
downwash  field  of  T2  is  exactly  equal  to  the  drag  of  r2  arising  from 
the  downwash  of  f1# 
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Consider  now  an  initial  distribution  designed  to  achieve  minimum 
drag.  (See  fig.  1.)  The  drag  added  by  a  small  additional  loading  will 
be  composed  of  three  parts,  namely: 

1.  The  drag  of  the  additional  lift  acting  alone 

2.  The  drag  of  the  original  loading  arising  from  the  downwash 

field  of  the  additional  loading 

3.  The  drag  of  the  additional  loading  induced  by  the  downwash  field 

of  the  original  loading 


Item  1  is  of  second  order  in  terms  of  the  magnitude  of  the  added  lift  for 
smooth  distributions,  that  is,  so-called  "weak  variations."  (The  fact 
that  this  second— order  term  is  invariably  positive  insures  that  the  drag 
will  be  a  minimum  and  not  a  maximum. )  Items  2  and  3  are  equal  by  the 
mutual  drag  theorem.  The  first— order  variation  in  drag  can  then  be  com¬ 
puted  by  considering  only  the  drag  of  the  small  additional  lift  acting  in 
the  induced  downwash  field  w-[(y)  of  the  original  lift. 

The  conditions  of  fixed  bending  moment  and  fixed  total  lift  are  met 
by  allowing  only  those  curves  of  lift  variation  that  produce  no  change  in 
these  quantities,  that  is,  curves  having  zero  area  and  zero  moment.  It 
can  be  seen  that  such  curves  of  variation  must  have  at  least  three  ele¬ 
ments  to  meet  the  conditions  of  zero  area  and  zero  moment.  Furthermore, 
any  curve  meeting  these  conditions  can  be  subdivided  into  groups  of  three 
elements  so  that  the  individual  groups  also  satisfy  the  conditions. 

Hence,  as  the  representative  of  such  restricted  curves  of  variation  we 
may  adopt  three  small  elements  having  areas  Zi,  Z2j  and  Z3  (fig.  l). 
These  elements,  together  with  their  positions  y^,  y2,  and  y3  and  the 
local  values  of  the  downwash  Wx1,  etc.,  due  to  the  original  loading 
must  satisfy  the  following  three  equations: 


for 

o 

It 

60 

Zi  +  =  0 

for 

o 

II 

PQ 

60 

li7i  +  ^2 72  +  l373  =  0 

(8) 

for 

SDi  =  0, 

liv±1  +  *2Wi2  +  l3v±  =  0 

j 


It  can  be  seen  that  these  equations  will  be  consistent  if  ~  a+byi, 

wi2  ~  a+by2  and  wi3  ~  a+by3,  where  a  and  b*  are  constants  to  be 
determined  from  the  given  conditions.  Since  such  equations  must  be  satis¬ 
fied  for  all  positions  yx,  y2,  etc.,  it  is  concluded  that,  in  general, 


wi  ^  a  +  by 


(9) 
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Hence,  for  a  minimum  induced  drag  with  a  given  total  lift  and  a  given 
bending  moment  the  downvash  must  show  a  linear  variation  along  the  span.2 
(See  fig.  2.) 

The  foregoing  method  may  be  readily  extended  to  a  more  general  class 
of  problems  involving  bending  moments  or  rolling  moments.  Suppose,  for 
example,  a  braced  wing  is  considered,  as  in  the  dotted  outline  of  fig¬ 
ure  3.  In  this  case  the  bending  moment  developed  by  that  portion  of  the 
lift  acting  inboard  of  the  point  of  bracing  attachment  may  be  of  no  con¬ 
cern,  but  it  may  be  desired  to  limit  the  bending  moment  developed  by  that 
portion  of  the  spanwise  load  curve  extending  between  this  point  and  the 
tip.  In  this  case  s0  will  not  be  zero.  At  least  three  elements  are' 
required  to  preserve  stationary  values  of  the  lift  and  bending  moment, 
and  it  is  evident  that  at  least  two  of  the  elements  must  lie  to  the  right 
of  the  point  sQ.  The  three  simultaneous  equations  are  (see  fig.  2): 

x  ^2  ^3  =  0 

^2(y2~s0)  +  23(y3— s0)  =  0 

liVi  +  ^2Wi2  +  ^3vi  =  0 

Here  y2  and  y3  are  to  the  right  of  the  point  sQ  and  j±  lies  to  the 
left  of  this  point.  For  these  equations3  to  be  consistent  wi  must  have 
the  form 


(10) 


wil  ~  aJ  wi2  ~  a+t>(y2~So)i  Wig  <v  a+b(y3-So) 


Hence ,  in  general,  the  downvash  will  he  a  constant  over  the  portion  of  the 
span  for  which  the  moment  is  not  specified,  as  illustrated  in  figure  3. 

If  no  restriction  whatever  is  placed  on  the  moment  there  is  obtained  the 
solution  of  Muni's  original  problem,  namely,  that  the  downvash  should  be 
constant  over  the  entire  span. 


p  . . . 

It  may  be  noticed  at  this  point  that,  whereas  the  discussion  bag  enrpha— 
sized  the  idea  of  minimizing  the  drag,  the  analysis  actually  makes  no 
distinction  between  the  lift,  bending  moment,  or  drag,  in  that  station¬ 
ary  values  of  all  three  are  demanded.  Thus  equation  (9)  may  be  consid¬ 
ered  a  necessary  condition  for  the  solution  of  the  following  problems: 
(l)  given  the  total  lift  and  the  induced  drag  to  find  the  distribution 
of  lift  over  the  span  that  will  result  in  a  minimum  bending  moment, 

311,1  (2)  given  the  bending  moment  and  the  induced  drag  to  find  the  dis¬ 
tribution  resulting  in  the  mazimum  total  lift. 

3See  reference  3  for  a  discussion  of  solutions  of  such  equations. 
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Determination  of  Span  Loading  and  Induced 
Drag  From  the  Downwash  Distribution 


The  case  of  bilateral  symmetry  with  moment  specified  about  the 
root  section  will  serve  as  an  example  of  the  calculation  of  the  actual 
span  loading  and  induced  drag.  It  will  be  evident  from  the  foregoing 
that  the  downwash  distribution  will  consist  of  two  straight-line 
segments  with  a  reversal  of  slope  at  the  plane  of  symmetry.  It  is  then 
necessary  to  compute  the  spanwise  variation  of  T  corresponding  to 
such  a  curve  of  downwash. 

To  perform  this  calculation  by  standard  methods  of  airfoil  theory, 
use  is  made  of  the  idea  that  at  a  great  distance  behind  the  wing  the 
vortex  sheet  forms  a  two-dimensional  field  of  motion,  with  the  discord- 
tinuity  in  the  lateral  velocity  across  the  sheet  given  by  dr/dy,  and 
the  downwash  w  given  by  twice  the  value  of  the  induced  downwash  wj 
at  the  wing.  Hence,  the  quantity  l/2  (dr/dy)  —  2iw^  can  be  evaluated 
by  means  of  the  familiar  complex  velocity  function  v  —  iw  of  the  two- 
dimensional  potential  theory  using  for  v  its  value  just  above  the 
vortex  sheet.  In  this  theory  if  the  vertical  component  of  velocity  w 
is  given  along  the  line  representing  the  trace  of  the  span,  then  the 
velocity  vector  at  any  other  point  in  the  field  £  =  y  +  iz  may  be 
obtained  from  the  relation  (reference  4) 


v  —  iw  = 


1  1  p+  a 

V  S2  -  P 


w(t)) 

tF| 


dri 


As  noted  above. 


dr 

5y 


=  v(y  +  oi)  -  v(y 


=  2v(y  +  oi) 


(ID 


(12) 


Introducing  v  =  a  +  by  for  y  >  0  and  v  =  a  -  by  for  y  <  0  into 
equation  (ll)  yields,  after  integration, 


i 

» 


and  hence 


—  =  -2a 


-=£=  +  -*> 
</  s2— y2  « 


cosh" 

s 


-JL- 


(13) 


A  2b  o  , -l 

+  —  y^  cosh 

JT 


The  spanwise  loading  thus  contains  the  elliptical  distribution  as  one 
component . 
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Equation  (13)  for  the  spanwise  distribution  of  circulation  enables 
the  determination  of  the  over— all  lift,  bending  moment,  and  drag  in  terms 
of  the  unassigned  constants  a  and  b.  The  use  of  equations  (2),  (6), 
and  (7),  together  with  the  wing  semispan  s,  yields  the  following  values: 


L  =  pVs‘ 


ita  +  bs 


B  =  pVs3  (  |  a  +  -  bs 

3  * 


(14) 


It  is  convenient  to  specify  the  bending  moment  of  the  lift  in  terms  of 
the  lateral  position  of  the  centroid,  or  center  of  pressure,  of  the  load 
curve.  The  lateral  centroid  as  a  fraction  of  the  semispan  s  may  be 
denoted  by  y*  (i.e.,  y!  =  2B/Ls).  Then,  solving  for  a  and  b. 


a 


> 


(15) 


The  expression  for  induced  drag  in  terms  of  the  lift  and  the  lateral 
center  of  pressure  becomes 


(2s) 


* 


12  ny* 


(16) 


This  equation  yields  the  minimum  drag  for  the  given  position  of  y1  .  If 
the  lateral  center  of  pressure  is  specified  so  as  to  coincide  with  that 
for  an  elliptical  loading  (i.e.,  b  =  0;  y'  =  4/3*),  then  the  above 
formula  reduces  to 

D±  =  — — -  (17) 

*|v2(2s)2 


The  optimum  distribution  of  loading  for  a  given  position  of  the  centroid 
y’  may  be  obtained  from  equation  (13)  with  the  aid  of  equations'  (15). 
The  result  is 
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Drag  For  a  Given  Bending  Moment  With  Unrestricted  Span 


The  foregoing  calculations  show,  as  was  to  he  expected,  that  the 
elliptic  loading  yields  a  smaller  drag  than  any  of  the  others  within  a 
restricted  span.  However,  if  the  restriction  on  the  span  is  removed, 
still  lower  values  of  the  induced  drag  can  he  obtained  without  any 
increase  in  the  bending  moment  at  the  wing  root.  The  lower  values  are 
obtained  by  permitting  the  span  to  increase  and  at  the  same  time  adopt¬ 
ing  a  more  tapered  form  of  the  loading  curve. 

Equation  (l6)  which  contains  the  three  variables  lift,  span,  and 
center  of  pressure  can  be  easily  rearranged  to  show  the  variation  of 
drag  with  span  when  the  bending  moment  and  the  lift  are  held  at  fixed 
values.  In  this  case,  the  lateral  position  of  the  center  of  pressure 
y1  s  will  be  fixed,  while  the  form  and  extent  s  of  the  load  curve  will 
vary.  In  order  to  provide  a  convenient  basis  for  comparison  the  span 
and  shape  of  the  load  curves  will  be  related  to  the  elliptic  loading. 

If  s/se  denotes  the  ratio  of  the  semispan  of  the  wing  to  that  of  an 
elliptically  loaded  wing  having  the  same  total  lift  and  bending  moment,' 
then  equation  (16)  can  be  rewritten: 


D,  = 


*  £  V2(2se)2  L 


Sp\4 

8  (  -2. 


16  (t)°  *  9 


(19) 


The  quantity  in  the  bracket  is  the  ratio  of  the  induced  drag  to  that  of 
the  corresponding  elliptically  loaded  wing.  This  ratio  is  plotted  in 
figure  4  to  show  the  decrease  of  drag  possible  by  increase  of  the  span. 
The  forms  of  load  curve  required  for  the  minimum  drag  at  various  values 
of  s/se  are  shown  in  figure  5- 


It  will  be  noted  that  a  15— percent  reduction  of  the  induced  drag 
below  that  for  elliptic  loading  can  be  achieved  with  a  15— percent 
increase  in  span.  Further  increases  of  span  between  15  percent  and 
50  percent  (s=1.15  to  1.50)  yield  no  significant  reductions,  however. 

At  still  larger  values  of  s  the  drag  becomes  lower,  and  approaches  zero 
at  an  infinite  value  of  s.  For  extreme  values  of  s/se  the  curves  begin 
to  show  negative  loadings  at  the  tips  and  eventually  the  bending  moment 
at  certain  points  along  the  span  will  exceed  that  at  the  wing  root. 


Ames  Aeronautical  Laboratory, 

National  Advisory  Committee  for  Aeronautics, 
Moffett  Field,  Calif.,  Sept.  25,  1950. 
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APPENDIX 

DEFINITIONS  OF  SYMBOLS 

L\  total  lift 
l  element  of  lift 

D^  induced  drag 

B  bending  moment 

p  air  density 

T  circulation 

induced  downwash  velocity  at  wing 
w  downwash  velocity,  at  infinity  (  w  =  2w±) 
v  lateral  velocity 

V  velocity  of  flight 

y,T|  distances  along  wing  semi  span 
sQ  point  of  origin  for  bending  moment 

s  length  of  wing  semi  span 

y*  lateral  position  of  load  centroid  as  a  fraction  of 


a,b  constants 
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FIGURE  I-  SPANWISE  LOAD  CURVE  WITH  THREE 
ELEMENTS  OF  VARIATION 


FIGURE  2r  DOWNWASH  DISTRIBUTION  FOR 
MINIMUM  DRAG  WITH  RESTRICTED  BENDING 
MOMENT  ABOUT  CENTER  SECTION 
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FIGURE  4.-  VARIATION  OF  INDUCED  DRAG  WITH  SEMISPAN  RATIO 
S/Se»  FOR  WINGS  HAVING  A  FIXED  TOTAL  LIFT  AND  A 
FIXED  BENDING  MOMENT. 
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NACA-Langley  -  12-21-50  -  1075 


FIGURE  5.-  VARIATION  OF  SHAPE  OF  THE  SPANWISE  LOADING 
CURVE  WITH  SEMISPAN  RATIO  S/Se,  FOR  WINGS  HAVING 
A  FIXED  TOTAL  LIFT  AND  A  FIXED  BENDING  MOMENT 
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